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ABSTRACT 

We analyze the processes relevant for star formation in a model with dark 
matter in the form of sterile neutrinos. Sterile neutrino decays produce an x- 
ray background radiation that has a two-fold effect on the collapsing clouds of 
hydrogen. First, the x-rays ionize the gas and cause an increase in the fraction 
of molecular hydrogen, which makes it easier for the gas to cool and to form 
stars. Second, the same x-rays deposit a certain amount of heat, which could, 
in principle, thwart the cooling of gas. We find that, in all the cases we have 
examined, the overall effect of sterile dark matter is to facilitate the cooling of 
gas. Hence, we conclude that dark matter in the form of sterile neutrinos can 
help the early collapse of gas clouds and the subsequent star formation. 

Subject headings: dark matter, neutrinos, reionization, star formation 

1. Introduction 

Numerous observations, from galaxy rotation curves to gravitational lensing, to cosmic 
microwave background radiation, all point to the existence of dark matter which is not made 
of ordinary atoms, but, rather, of some new, yet undiscovered particles. In particular, the 
abundance of dark matter has been precisely determined recently by the WMAP (Spergel 
et al. 2003, 2006). However, the nature of dark-matter particles remains unknown. One 
attractive candidate is a sterile neutrino with mass of several keV and a small mixing with the 
ordinary neutrinos. If such a particle exists, it could be produced in the early universe with 
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the right abundance to be dark matter either from neutrino oscillations (Dodelson & Widrow 
1994; Shi & Fuller 1999; Abazajian, Fuller & Patel 2001; Abazajian, Fuller & Tucker 2001; 
Abazajian 2006a; Dolgov & Hansen 2002), or from some other processes, for example from 
inflaton decay (Shaposhnikov & Tkachev 2006). The same particle can explain the observed 
velocities of pulsars because its emission from a cooling neutron star would be anisotropic, 
hence providing the neutron star with a recoil velocity of a sufficient magnitude (Kusenko 
& Segre 1997, 1999; Kusenko, Segre & Vilenkin 1998; Fuller et al. 2003; Barkovich, D'Olivo 
& Montemayor 2004, 2005; Kusenko 2004). In addition, the neutrino kick can enhance 
the convection during the first second of the supernova explosion, which can increase the 
energy of the supernova shock and bring the supernova calculations in better agreement with 
the observations (Fryer & Kusenko 2006). Sterile neutrinos can also help the formation of 
supermassive black holes in the early universe (Munyaneza & Biermann 2005; Fan et al. 
2001). 

Although the dark-matter sterile neutrinos have interactions that are too weak to be 
discovered in the laboratory, some other sterile neutrinos may exist. Most models of neu- 
trino masses introduce sterile (or right-handed) neutrinos to generate masses of the ordinary 
neutrinos via the seesaw mechanism (Gell-Mann et al. 1979; Yanagida 1980; Mohapatra & 
Senjanovic 1980). The name sterile was coined by Pontecorvo (1967). Many seesaw models 
assume that sterile neutrinos have very large masses, which makes them unobservable. How- 
ever, a number of theoretical models predict the existence of sterile neutrinos at or below 
the electroweak scale (Farzan et al. 2001; Appelquist & Shrock 2002; Asaka & Shaposhnikov 
2005; de Gouvea 2005). There is an indication that some fight sterile mass eigenstates may 
cause the neutrino flavor transformation seen in the Los Alamos Liquid Scintillator Neutrino 
Detector (LSND) experiment (Athanassopoulos et al. 1998; Sorcl, Conrad & Shaevitz 2004; 
dc Gouvea 2005). The results of LSND are being tested by the MiniBooNE experiment. 
One can expect that there arc at least three light sterile neutrinos (Asaka, Blanchet & Sha- 
poshnikov 2005), in which case the neutrino oscillations can explain the baryon asymmetry 
of the universe (Akhmedov, Rubakov & Smirnov 1998; Asaka & Shaposhnikov 2005). 

Although dark-matter sterile neutrinos are stable on cosmological time scales, they nev- 
ertheless decay. The mass matrix is assumed to have the seesaw form with the Majorana 
mass term of the order of several keV and with a much smaller Dirac mass term. The diag- 
onalization of such a mass matrix yields the light mostly active Majorana mass eigenstates 
and the heavy sterile Majorana eigenstates, which can decay into the light ones. The dom- 
inant decay mode, into three light neutrinos, is "invisible" because the daughter neutrinos 
are beyond the detection capabilities of today's experiments. The most prominent "visible" 
mode is decay into one active neutrino and one photon, i/g — > i^al- Assuming two-neutrino 
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mixing for simplicity, one can express the inverse width of such a decay as 
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where is the mass and 6 is the mixing angle. For mass units we set c = 1. 

Since this is a two-body decay, the photon energy is half the mass of the sterile neutrino. 
The monochromatic line from dark matter decays can, in principle, be observed by x-ray 
telescopes. No such observation has been reported, and some important limits have been 
derived on the allowed masses and mixing angles (Abazajian, Fuller & Patel 2001; Abazajian, 
Puller & Tucker 2001; Boyarsky et al. 2005, 2006a,b,c; Riemer-S0rensen et al. 2006; Abazajian 
2006b; Watson et al. 2006). These constraints are based on different astrophysical objects, 
from Virgo and Coma clusters, to Large Magellanic Clouds, to the Milky Way halo and its 
components, etc. There are different uncertainties in modeling the dark matter populations 
in these objects. Different groups have also used very different methodologies in deriving 
these bounds: from a conservative assumption that the dark-matter line should not exceed 
the signal, to more ambitious approaches that involved modeling the signal or merely fitting 
it with a smooth curve and requiring that the line-shaped addition not affect the quality of 
the fit. The x-ray fimits are usually shown based on the assumption that the sterile neutrinos 
constitute the entire dark matter (fig = ^dm, where s and dm refer to the sterile neutrinos 
and the dark matter, respectively). However, the x-ray observations exclude certain masses 
and mixing angles even if the sterile neutrinos make up only a fraction of dark matter. The 
production via mixing (Dodelson & Widrow 1994) provides the lowest possible abundance 
of sterile neutrinos. In Fig. 1 we show both the bounds based on the assumption Qg = Qa-m 
and the model-independent exclusion region (Kusenko 2006) based on the production only 
via the Dodelson & Widrow (1994) mechanism. Kusenko (2006) used the analytical fit to 
the numerical calculation of Abazajian (2006a), and re-scaled the x-ray fiux in proportion 
to the amount of relic neutrinos produced for a given mass and mixing angle. 

A different set of bounds comes from the observations of the Lyman-alpha forest (Viel 
et al. 2005, 2006; Seljak et al. 2005), which limit the sterile neutrino mass from below. Based 
on the high-redshift data from SDSS and some modeling of gas dynamics, one can set a limit 
as strong as 14 keV (Seljak et al. 2005). However, the high-redshift data have systematic 
errors that are not well understood (Jena et al. 2005), and more conservative approaches, 
based on the relatively low-redshift data, have led to some less stringent bounds (Viel et al. 
2005). More recently, Viel et al. (2006) have reanalyzed the high-redshift data and arrived at 
the bound > 10 keV. The mass bound as quoted depends on the production mechanism 
in the early universe. 



The Lyman-alpha observations constrain the free-streaming lengths of dark matter par- 
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ticles, not their masses. For each cosmological production mechanism, the relation between 
the free-streaming length and the mass is different (Hansen et al. 2002). For example, the 
bound > 10 keV (Viel et al. 2006) applies to the production model due to Dodelson & 
Widrow (1994). Even within a given cosmological scenario, there are uncertainties in the 
production rates of neutrinos for any given mass and mixing angle (Asaka, Laine & Shaposh- 
nikov 2006). These uncertainties may further affect the interpretation of the Lyman-alpha 
bounds in terms of the sterile neutrino mass. The x-ray bounds depend on a combination of 
mass and mixing angle as in eq. (1), regardless of one's cosmological assumptions, as long as 
all the dark matter is comprised of sterile neutrinos. However, the x-ray bounds are subject 
to uncertainties in modeling the dark matter populations of different objects. 

If the sterile neutrinos make up only a part of dark matter, none of the above bounds 
apply. Also, if inflation ended with a low reheating temperature, the bounds are significantly 
weaker (Gelmini, Palomares-Ruiz & Pascoli 2004). 

It should also be mentioned that the Lyman-alpha bounds appear to contradict the 
observations of dwarf spheroidal galaxies (Wilkinson et al. 2006; Gilmore 2006; Strigari 
et al. 2006), which suggest that dark matter is warm and which would favor the 1-5 keV 
mass range for sterile neutrinos. There are several inconsistencies between the predictions of 
N-body simulations of cold dark matter (CDM) and the observations (Willman et al. 2004; 
Bode, Ostrikcr & Turok 2001; Peebles 2001; Dalcanton & Hogan 2001; Bosch & Swatcrs 2001; 
Zentncr & Bullock 2003; Abazajian, Koushiappas & Zcntncr 2006; Sommer-Larscn & Dolgov 
2001; Govcrnato et al. 2004; Kormcndy & Fisher 2005; Fcllhaucr ct al. 2006; Belokurov et al. 
2006; AUgood et al. 2006). Each of these problems may find a separate independent solution. 
Perhaps, a better understanding of GDM on small scales will resolve these discrepancies. It 
is true, however, that warm dark matter in the form of sterile neutrinos is free from all these 
small-scale problems altogether, while on large scales WDM fits the data as well as GDM. 

Several independent considerations, from the pulsar kicks to small-scale behavior of 
dark matter, favor sterile neutrinos as the dark matter candidate. It is important, therefore, 
to examine the effects of sterile-neutrino dark matter on star formation and reionization 
of the universe. Yoshida et al. (2003) have attempted to study such effects for some kind 
of warm dark matter, but they have used a linear matter power spectrum and assumed 
that the dark matter particles have a thermal distribution, which is not the case for sterile 
neutrinos (Abazajian 2006a). Even more importantly, they have not considered the effects 
of sterile neutrino decays and the resulting x-ray background. It was shown by Biermann & 
Kusenko (2006) that the x-rays produced by the decays of relic sterile neutrinos can speed 
up the formation of molecular hydrogen, which can help the cooling of the gas and the star 
formation, which can lead, in turn, to a reionization of the universe at the redshift consistent 
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with the WMAP results. 

More specifically, the WMAP (three years) measurement of the reionization redshift 
Zr — lO.Q^Ig (Spergel et al. 2006) has posed a challenge to theories of star formation. On 
the one hand, stars have to form early enough to reionize gas at redshift 11. On the other 
hand, the spectra of bright distant quasars imply that reionization must be completed by 
redshift 6. Stars form in clouds of hydrogen which collapse at different times, depending on 
their sizes: the small clouds collapse first, while the large ones collapse last. If the big clouds 
must collapse by redshift 6, then the small halos must undergo the collapse much earlier. It 
appears that the star formation in these small halos would have occurred at high redshift, 
when the gas density was very high, and it would have resulted in an excessive Thompson 
optical depth (Haiman & Bryan 2006). To be consistent with WMAP, the efficiency for the 
production of ionizing photons in minihalos must have been at least an order of magnitude 
lower than expected (Haiman & Bryan 2006). One solution is to suppress the star formation 
rate in small halos by some dynamical feedback mechanism. The suppression required is by 
at least an order of magnitude. 

Alternatively, one can consider warm dark matter, in which case the number of small 
clouds is strongly suppressed. However, for some WDM candidates (for example, the grav- 
itino), the lack of power on small scales in the spectrum of density perturbations leads to an 
unacceptable delay in the onset of star formation (Yoshida et al. 2003). This is not the case 
for sterile neutrinos, which decay slowly during the " dark ages" preceding the birth of stars. 
The decay photons could ionize gas, and, although this ionization is too weak to affect the 
WMAP measurements directly (Mapelli & Ferrara 2005; Mapelli et al. 2006), the increase in 
the ionization fraction could cause a rapid production of molecular hydrogen at redshifts as 
high as 100 (Biermann & Kusenko 2006). Molecular hydrogen is the essential cooling agent 
that allows the gas to cool and collapse to form the first stars. 

Biermann & Kusenko (2006) have not considered the effects of gas heating associated 
with the decays, which may, in some cases, thwart the effects of ionization. The increase in 
the fraction of molecular hydrogen can help coohng the gas, but if the heating from sterile 
neutrino decays is significant, it can overcome the radiative cooling and stymie the collapse 
and star formation. In this paper we will examine thermal evolution of gas cloTids in detail, 
taking into account both effects of the sterile neutrino decays, namely the ionization and the 
heating of gas. 

The paper is organized as follows. In Section 2 we will discuss the role of sterile neutrino 
decays and the production of molecular hydrogen in the cooling of the primordial gas clouds. 
Next, in Sections 3-6 we will brieffy describe the top-hat overdensity simulation. The results 
are shown in Section 7. 
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2. H2 cooling and sterile neutrinos decays 

The existence of the first objects is a direct consequence of the growth of the primor- 
dial density fiuctuations. At the beginning, there are hnear density perturbations which 
expand with the overall Hubble flow. Subsequently, these perturbations can grow and form 
primordial clouds. Clouds with enough density contrast decouple from this flow and start to 
collapse. The kinetic energy of the infaUing gas is dissipated through shocks and the cloud 
becomes pressure supported. The subsequent evolution of the cloud is determined by its abil- 
ity to cool sufficiently fast. Clouds that do not cool fast enough stay in a pressTire-supported 
state and do not form stars. The existence of the efficient cooling mechanism is necessary 
to continue the collapse of the cloud, its subsequent fragmentation and star formation. 

In the absence of metals, the most important cooling mechanism is cooling by hydrogen 
atoms and molecules. This mechanism involves the collisional excitation of H or Ha, sub- 
sequent spontaneous de-excitation, and photon emission. The emitted photons can escape 
and take away energy, hence reducing the kinetic energy of the collapsing cloud. Cooling 
by H dominates for temperatures above 10^ K, while the H2 cooling is more important at 
lower temperatures. The temperature of the primordial cloud can be higher than 10"^ K only 
during a short period of virialization of the baryonic gas. Therefore, H2 cooling is the most 
important cooling mechanism during almost entire evolution of the halos, and it is crucial 
for the star formation. 

In primordial gas clouds, H2 molecules can form mainly through the coupled reactions 

e- + H H-+'y (2) 
H- + H ^ H2 + e-, (3) 

or 

H+ + H H+ + J (4) 
H+ + H H2 + H+, (5) 

in which e~ and H'^ act only as a catalyst. X-ray radiation can increase the production of 
the H2 by enhancement of the ionization fraction. 

As discussed above, sterile neutrinos with the mass of several keV and a small mixing 
are stable on cosmological time scales. Nevertheless some of them can decay. The decay 
channel important for us is 

'^s^l^a + 7- (6) 

Its inverse width is given in equation (1). The energy of the photon is equal to 

_ (mg + ml) rris 
- T' 
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where rua is the mass of the active neutrino. X-ray photons from the sterile neutrinos decays 
can increase the ionization fraction and the H2 production. This can subsequently lead to a 
speed-up of the gas cooling and star formation. 



3. Density evolution 

We will follow the evolution of the baryonic top-hat overdensity, the gas temperature 
and its H2 and e" fraction. Our goal is to juxtapose the evolution of the gas temperature in 
the primordial clouds in the ACDM model and the WDM model with keV sterile neutrinos. 
We assume that the dark matter is smoothly distributed in the region we are considering. 

We will use the parameters from the best fit to the three years WMAP data (Spergel et 
al. 2006), namely Ho = 100/i km s'^ Mpc'^, h = 0.73, Qq = ^dm + ^b^ 0.238, = 0.042, 
^dm ^Qo-^b^ 0.196, Qa = 0.762, Q^* = Qam + ^6 + ^^a = 1- 

Let us take some spherical region with the top-hat overdensity 5. It is well known that 
such an overdensity will evolve according to the following equation (Padmanabhan 1993) 

l + ^^9 (a-sinc.)^ ^ (8) 
2 (1 — cos a) 

where the parameter a is related to the redshift z and redshift of virialization Zyir through 

2/3 

(9) 



1 + Zr ir / a — sin a 



27r 

Further evolution depends on the type of the matter in the overdense region. If it is 
the dark matter then after virialization its density remains constant forever. The situation 
is different in the case of baryons. If there exists an efficient cooling mechanism, then the 
density gradually increases. Otherwise, the density remains constant, and there is no star 
formation in the halo. 

In reality an overdense region consists of two components namely, baryonic gas and dark 
matter. Baryonic density fluctuations are suppressed before recombination. In contrast, the 
dark matter density fluctuations can start to grow well before this period, so its initial 
overdensity is larger than that of the baryons. Dark matter forms the potential wells into 
which the baryonic matter falls. Overall, the baryon density follows the dark matter density 
in the early universe. Therefore, one can consider the primordial fluctuations as the halos 
consisting of baryonic and dark matter, whose overall overdensity follow the top-hat solution 
(eqs. 8 and 9). Of course, this assumption breaks down at the later stages of the halo 
evolution, e.g. after the virialization of the dark matter. 
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FoUowing Tegmark et al. (1997), let us assume that the evolution of the baryonic and 
dark matter density is described by equations (8) and (9) until a reaches 37r/2. This occurs 
at the redshift 

= 1-06555 (1 + z,ir) - 1- (10) 

After this redshift, we assume that the density of the halo Pb+dm stays constant e.g. it is 
equal to the value 

Pb+dm = ISTT^f^oPo (1 + Z^irf , (H) 

where po is the critical density of the universe. 

We note that the overdensity 5 in equation (8) is given by 

(5=^^-1, (12) 
P 

where p is the combined average density of dark and baryonic matter. 

We have taken the ratio of the baryonic to the dark matter mass in the halo to be the 
standard cosmological ratio of O^/^dm ~ 0.214. The density of the baryonic matter can be 
expressed as 

Pb — Pb+dmT^- (13) 
I In 



4. Temperature evolution 

The gas temperature evolution is governed by the equation of energy conservation 

du _ p_dpb _ A_ 
dt pI dt pb 

where p, u and A are the pressure, the internal energy per unit mass and the cooling/heating 
function respectively. We use the equation of state of perfect gas 

P=(7-1)P6«, (15) 

where 7 is the adiabatic index. 

The internal energy per unit mass can be described as 

(16) 

7-1 pmH 
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where k is the Boltzmann constant, T is the gas temperature, denotes the molecular 
weight, and uih is the mass of the hydrogen atom. Introducing the number density of non 
dark matter particles Up we can write pb as 

Pb = pmHUp. (17) 

If we use eqs. (15), (16) and (17) than we can cast eq. (14) in the form 

dT , ( T dn,o A \ T da T dj 

— = (7 - 1) £ +7 c _| L_ (18) 

dt \np dt Upk J p dt (7 — 1) dt 



The time t can be translated to the redshift z in the following way 

dt _ 1 

Ho (1 + z) ^JnA + no{l + zf 

Using eqs. (18) and (19), we obtain 



(19) 



dT , . T dur, T da T d'y (7 — 1) A . . 

d; = ^^-^^^^^^7^Tz^h^)Tz^ / (20) 

To calculate the gas temperature of the primordial cloud we have to integrate this equation. 

One can assume that the gas temperature evolution is governed by eq. (20) except for 
the redshift range between Z2t^/2 and z^ir- If we simply integrate eq.(20) to the redshift z^r, 
than the gas temperature will be much lower than the virial temperature T^ir, at least for 
clouds with mass greater than some critical value (Kitayama et al. 2001) 

r.. = 9.09 X 10' K (^) (^i|5^) (%r)) (1 + , (21) 

where p^r, Mb+dm a-nd Ac {z^ij.) are respectively the mean molecular weight during virial- 
ization, the halo mass and the overdensity of the halo at virialization given by (Bryan & 
Norman 1998), which is very well approximated by ISvr^ for redshifts z^r < 100. Therefore 
we have to take into account shocks and increase the gas temperature to the virial value. 
It is reasonable to assume that the gas temperature evolution is linear between the redshift 
z^T^l2 and Zyir and is described by the equation 



dT _ T^i-f. — T37r/2 
dz Zyij. — ^37r/2 

where T3^/2 is the temperature at the redshift 2;37r/2- 



(22) 
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Let us consider the following five species: H, H+, H2, e , and He with the mass fraction 
Y = 0.244 (Izotov & Thuan 1998). Let us denote by (i =H, H+, H2, e" and He) the 
number density of baryonic component i. We assume that nH+ = Ug e.g. there is no He or 
H2 ionization. The number density of hydrogen can be written as 

uh + nH+ + '2nH2 ■ (23) 
If we introduce the number fraction of component "i" 

Xi = ^, (24) 
n 

then the particle number density Hp and the baryon density ph are equal to 

np = n(^l + ^-XH2+ (25) 

and 

Pb = -j^, (26) 
where X = l — y = 0.756 is the hydrogen mass fraction. 
Using eqs. (17), (25) and (26) we can get 

At the virialization we have typically ~ and ~ 1, so one can take the mean molecular 
weight at this time to be Hvir = 1/ (2-^^ + 0.25y) 0.63. We can express 7 as 



where 77? = 7^?+ = 7e = ^He = 5/3 and -fH2 = 7/5. 



5. Chemical reactions and cooling/heating function 

To calculate the coohng/heating function A in eq. (20) we have to take into account 
all of the relevant chemical and thermal processes and follow the number density evolution 
of different baryonic components. All chemical and thermal processes included in our cal- 
culation are summarized in Tables 1 and 2. The evolution of Xg and Xh+ is described by 
equations 

dx 

—^^{A (z) + ki2 + kuUe) xh - haxln (29) 



- 11 - 



— ^ = 30) 

dt dt' ^ ' 

where A {z) is the production rate of free electrons by the sterile neutrino decays which is 
calculated in the appendix. In eq. (29) Only the ionization and recombination of hydrogen 
are included. All the rates ki are given in Table 1. 

We note in passing that neglecting He and H2 ionization renders the system with fewer 
free electrons, so the production of H2 is slower. In addition, cooling by is stronger 
than by H2 (see last figure in Galli & Palla (1998)). Therefore, the cooling we get in this 
approximation is somewhat less efficient than it should be in reality. 

The molecular hydrogen number fraction can be calculated by integration of the follow- 
ing equation 



dXH2 

dt 



where 



= Uh {Xekm - k-iXH^) " {hXH2 + h^e) riHi " (^^10 + hi) Xh2, (31) 

, kik2 , k4k^ . . 

km = 7 — —j—/ — + -r — —j—, — 32) 

k2 + k3/nH k^ + ke/nn 

is the molecular hydrogen formation rate via H~ and channel (Tegmark et al. 1997). 
To get eq. (32) we have to make the assumption that after formation of H" and H2 they 
are used almost instantaneously in the next reactions (reactions number 2 and 3 or 5 and 6 
respectively, see Table 1). This is a very good approximation, because in reahty the H~ and 
H2 fractions always remain small. 

The cooling/heating function A is given in the appendix (eq. C3). It includes contribu- 
tions from the following processes: 



1. excitation of H by collisions with e~ 

2. recombination to H 

3. ionization of H by collisions with e~ 

4. photo-ionization of H by the CMB radiation 

5. excitation of H2 by collisions with H and H2 

6. H2 formation cooling 

7. dissociation of H2 by collisions with H, H2 and e 

8. photo-dissociation of H2 by CMB photons 
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9. Compton cooling 

10. bremsstrahlung 

11. photo-ionization of H by the photons from the sterile neutrino decays and H ionization 
by the secondary electrons produced due to these photons 

We have not inchidcd the effect of the photo-dissociation of H2 by the photons from the 
sterile neutrino decays. This is because the flux of these photons is small in comparison 
with the CMB radiation in the 10 - 30 eV range, which is the energy range relevant to the 
photo-dissociation of H2. 

The formulae for the contributions to the A function from all of these processes are 
given in Table 2. We have assumed the ortho-H2 to para-H2 ratio to be 3 : 1. 

6. Initial conditions and numerical procedure 

We start our calculations at the redshift of recombination Zrec = 1013, taking as initial 
conditions Xg = 1, = and assuming that the gas temperature is equal to the CMB 
temperature. Evolution of the halo is insensitive for the initial value of xh2- 

We follow the evolution of the density, eqs. (8) and (9), and the temperature, eq. (20) 
with the cooling/heating function given in eq. (C3) untill the redshift is equal to 2;37r/2, 
cf. eq. (10). From that time we assume that the density of the halo stays constant and is 
equal to its virial value (11). We assume that the gas starts to be shock heated in this period 
and reaches its virial temperature at z^r. The temperature of the gas is assumed to change 
linearly with the redshift between Z3.„/2 and z^r, eq. (22). For the redshift lower than z^r 
the gas temperature is calculated once again by integration of eq. (20). 

Following Tegmark et al. (1997), we use the following criterion for the successful collapse: 

T {r]Zyir) < r]T (z^ir) , (33) 
where r] = 0.75 (Tegmark et al. 1997). 

7. Results 

The results of our analysis depend, of course, on the parameters we choose for the 
dark-matter sterile neutrinos. The x-ray flux, and the resulting increase in the ionization 



fraction depend on the combination (mf sin^ 9). The power spectrum on small scales, which 
determines the presence or absence of small clouds, depends on the free-streaming length, 
which is set by the mass (in every given cosmological scenario). 

Based in part on the existing limits, we have chosen four illustrative benchmark cases 
to illustrate the effects of ionization and heating due to sterile neutrino decays on the cloud 
collapse. 

We have compared the following four models: 

1. Cold dark matter with no additional radiation besides CMB (CDM). 

2. Warm dark matter (WDM) with the sterile neutrino mass = 25 keV and mixing 
angle sin^ ^ = 3 x lO'^^ (WDMl). 

3. WDM with = 15 keV and sin^ ^ = 3 x IQ-^^ (WDM2). 

4. WDM with rris = 3.3 keV and sin^ ^ = 3 x lO'^ (WDM3). 

The values of parameters chosen for these models are shown in Fig.l. 

Some comments are in order. The choice of WDMl is obvious because it satisfies all the 
constraints, including the x-ray bounds and the Lyman-alpha bounds. The masses as large 
as rris = 25 kcV cannot be constrained by cither Chandra or XMM, and the best constraint 
comes from HE AO, which is weaker because of its lower energy resolution (Boyarsky et al. 
2005, 2006a,b,c). Also, WDMl corresponds to a very small free-streaming length and, on 
large scales, it creates structure indistinguishable from CDM. The only difference between 
WDMl and CDM is the additional x-ray flux. Therefore, the comparison of CDM and 
WDMl cases provides a clean test of the effects of the increased ionization fraction (Biermann 
& Kusenko 2006). Also, WDMl provides a realistic set of parameters consistent with all 
the bounds, as well as the pulsar kick mechanism due to resonant active-sterile neutrino 
conversions (Kusenko & Segre 1997). 

WDM2 is also consistent with all the existing limits, including the most ambitious 
Lyman-alpha bounds (Seljak et al. 2005). (We remind the reader that Viel et al. (2006) have 
obtained a 50% weaker limit using the same data.) The x-ray background is lower in this 
case than in the case of WDMl. 

Finally, WDM3 corresponds to masses that would be ruled out by the "strong" Lyman- 
alpha bounds obtained by Seljak et al. (2005) and by Viel et al. (2006). However, one has 
to remember that the relation between the masses and the free-streaming length depends 
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on the production mechanism of sterile neutrinos. For example, the free-streaming length 
in the Shi & Fuller (1999) production mechanism is different from the Dodelson & Widrow 
(1994) mechanism. Even within the latter framework, additional entropy production can red- 
shift the sterile neutrino momenta and eliminate the Lyman-alpha bound (Asaka, Kusenko 
& Shaposhnikov 2006). If dark matter is produced from a coupling between the sterile neu- 
trinos and the infiaton (Shaposhnikov & Tkachev 2006), the free-streaming length is, again, 
different. WDM3 requires redshifting of the neutrino momenta by a factor ~ 3, as compared 
to the thermal distribution. This can be achieved, for example, by entropy production of 
the order of 10^, which is comparable to the effect of change in the number of degrees of 
freedom from the electroweak scale to zero temperature (Kusenko 2006). Hence, WDM3 
should not be considered an unrealistic scenario. In addition, WDM3 is of interest because 
it is consistent with the pulsar kick mechanism via off-resonant active-sterile neutrino oscil- 
lations (Fuller et al. 2003), even if this sterile neutrino is not the dominant component of 
dark matter. 

The results of the calculation are presented in Figures 2, 3 and 4. Several conclusions 
can be drawn from these results. We find that large clouds, with masses greater than 10^°Mq, 
collapse regardless of the details of cooling in all the models defined above. One could expect 
this behavior because the virial temperature of such a cloud is higher than lO"* K, so the 
H cooling mechanism is efficient enough and the additional H2 does not play an important 
role. 

There is, however, a dramatic difference when it comes to smaller clouds. For a virial- 
ization redshift z^r = 20, a cloud of 4 x lO^M© collapses in model WDMl, but not in the 
CDM or other WDM models. This illustrates the importance of the increased H2 fraction 
on the cooling of gas. 

In WDM2 and WDM3 model, because of the relatively low associated x-ray flux, the 
temperature evolution is practically unaffected by the photons from the dark matter decay. 

From Fig. 4 one can draw another conclusion. In the large clouds, all of the H2 molecules 
are destroyed during the virialization process. This is because of the very high virial tem- 
perature. In contrast, there is a drastic increase in H2 fraction during the virialization of the 
clouds with the mass of 4 x lO^M©. In both cases, most of molecular hydrogen is produced 
after virialization. 

For clouds with M = 4 x IO^Mq and z^r = 100 the critical II2 fraction (see eq. (11) 
of Tegmark et al. 1997) is reached at virialization. On the other hand, in the clouds having 
equivalent mass with Zvir = 20, the molecular hydrogen faction xh^ attains its virial critical 
value only for WDMl model. This happens Bit z — 7.5. 
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In all the WDM models the ionization and the molecular hydrogen fraction are enhanced 
in comparison with the CDM model, as suggested by Biermann & Kusenko (2006). In 
addition, we can see that in WDM models the temperature of the collapsing clouds soon 
after virialization is lower then that in the CDM case, which speeds up the collapse, at least 
for the low-mass clouds. 

We note in passing that the density of the cloud was kept constant after virialization, 
so the evolution of T, and is reliable for the redshifts close to virialization, but not 
for much smaller redshifts. 

It appears that WDM3 model is disfavored. The increase in the cooling rate is almost 
negligible in comparison with the CDM case, while the suppression on the small scales is 
very strong. Therefore, this leads to an unacceptable delay in the onset of star formation, 
as discussed by Yoshida et al. (2003). 

In the case of WDMl and WDM2 models, the cooling is enhanced by the additional 
fraction of molecular hydrogen. However, we cannot definitively prove that star formation 
occurs at the redshift consistent with WMAP because of the simplifying assumptions and 
limitations of our analysis. 

The main limitation is that the virialization redshift of a particular cloud is a fixed 
parameter in the simulation. In reality, it can be changed by the enhancement of the cooling, 
so the cloud can reach its virial stage much earlier. Moreover, the density of the cloud is 
kept constant after virialization, which is a simplifying assumption. We can not calculate 
the real redshift of the cloud collapse and onset of the star formation. As a result, we cannot 
calculate how many clouds collapse at a given redshift and we cannot calculate the redshift 
of reionization in WDM model rehably. Also, we have not included the additional x-rays 
from the decay of sterile neutrinos inside the collapsing clump. This radiation can cause 
an additional increase in ionization and heating, which could change our results. Also the 
effects of opacity can be important in large clouds, which we have not taken into account. 
These limitations will be addressed in the future work. 



8. Conclusions 

We have performed a detailed analysis of the cooling and collapse of primordial gas in 
the model with warm dark matter, taking into account both the increase in the fraction of 
molecular hydrogen (Biermann & Kusenko 2006) and the heating due to the sterile neutrino 
decay. As expected, the effect on the largest gas clouds is negligible, but the smaller clouds 
are, in fact affected. We have performed the analysis for some benchmark cases which arise 
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in realistic scenarios. For the largest clouds, the additional molecular hydrogen made no 
difference. For smaller clouds, especially for those with masses 10^ — 10^ Mq, the increase in 
the x-ray background made the successful collapse possible in cases where it could not occur 
in the absence of sterile neutrino decays. 
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Appendix 

A. Specific intensity of the photons from the sterile neutrino decays 

In this section we shall calculate the photon spectrum produced by the radiative decays 
of the sterile neutrinos (our derivation is similar to the one in Masso & Toldra (1999)). The 
dominant channel of radiative decays is through the following reaction 

T^s^T^a + l, (Al) 

where Vg is the sterile neutrino and Ua is an ordinary active neutrino. The photon energy is 
given by 

^ = K±!!!|) ^ rn^^ (A2) 



where and nia are respectively the mass of the sterile and an active neutrino. The 
cosmic expansion redshifts the photon energy. A photon that at redshift z has energy E was 
produced at redshift Zq given by 

l + zo = (l + ^)^. (A3) 

Let Fe {z) be the energy flux (in units of erg cm~^ s~^) at redshift z of photons with 
energy E produced by the sterile neutrino decays. The flux per unit energy and solid angle 
is given by 

d'FE jz) ^ d^F^jz) 
dEdQ dEdn ' ^ ^ 
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where F„ (z) is the photon flux at redshift z. It is related to the photon flux at emission 
redshift Zq through the foUowing equation 

where wc have inchidcd the factor of (^j-f^^ produced by the expansion of the universe. 
The photon density emitted at Zq is given by the usual decay law 

6nj (zq) = —Us (zq) , (A6) 

T 

where r is the mean sterile neutrino lifetime for the decay into an active neutrino and a 
photon (see eq. 1), Ug (zq) is the number density of the sterile neutrino at redshift Zq and St 
is given by 

c , dt , ^fo 1 dE 

^^•^ Ho (1 + ^o) V^A + (1 + ^o)' ^ ^ 



where H [zq] = Hoy + + Zq)^ is the Hubble expansion rate at photon emission 
redshift Zq. Writing everything together we obtain 

d'FEjz) _ 1 fl + zy n4zo)c 
dEdVL A-K \l + zo) tH{zo)' ^ ' 

Choosing the time tp <^t otherwise arbitrary we can write 

where t (zq) is the age of the universe at redshift Zo, tp = t (zp) and Ug (z) would be the 
number density of sterile neutrinos at redshift z if they did not decay. 



Assuming that all of the dark matter consists of sterile neutrinos, Ug (z) can be written 



as 



Us (z) = n° (1 + zf , (AlO) 

where 

n° = nam-^. (All) 
rus 

In the case of a flat universe flUed with matter and a nonvanishing cosmological constant 
the age of the universe at redshift z is given by (Masso & Toldra 1999) 



2 \/ni+ JnA + no{i + zf 

t (z) = — In y . (A12) 
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Taking tp — and substituting eq. (A9) into eq. (A8) we can get the specific intensity 
(in units of erg cm~^ s~^ sr~^ Hz~^) of the photons with frequency u from the sterile neutrino 
decays at redshift z 

We have to multiply this value by 6^"^"^=*^^'^"^ in order to correct specific intensity for the 
absorption by the background neutral hydrogen atoms. The optical depth is given by 

rabs{z,zo) = J ai2 {E {z')) UH {z') c—dz' , (A14) 

where ai2 {E {z')) and nn {z') are the cross section for the H ionization at energy E [z') and 
the hydrogen number density of the Universe at redshift z' respectively; E {z') is the energy 
of the photon at redshift z' , which was emitted in sterile neutrino decay at redshift Zq and 
is given by the equation 

E{z')^\^E,. (A15) 

i- + Zq 

The mean hydrogen number density of the Universe at redshift z' can be written as 

riH {z') = [1 - X, {z)] (1 + z'f , (A16) 
where n^/ is its present value and is equal to 

The cross section for the H ionization is given by 

ai2(£;(^')) = ^o(|^)', (A18) 

where ctq = 7.909 x lO'^^ cm'^ and hi^g = 13.6 eV. 

Neglecting the fl\ term in eq. (19), which is a very good approximation at high redshifts, 
and using eqs. (A3), (A15), (A16) and (A18) we can cast eq. (A14) into the form 

Tabs {z, zo) = ao (^^) ' (1 + zf nlcH^'n-"' £ [1 - {z')] (1 + z'y'^' dz'. (A19) 
As Xe (z), let us take 

a;, (^) = I ° ^ < ^'•^^ (A20) 

I i Z ^ Zfrer , 
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where Zrec is the recombination redshift. It gives us 

2 / u.u \ 3 



Tabs {Z, Zq) = -C7o (1 + zf^ U^cHq ^Oq 



1/2 



1 + ^ 



1 + min (2:0, Zrec) 



3/2 



(A21) 



Photons emitted before recombination are Compton scattered and thermahzed. For 
simphcity we assume that these photons are not visible. The final specific intensity is equal 
to 



hcns{z)e — e—'aba(^'''o) 







hv > 

hv < Eoj^ 



where 



t{zo)=t(^^{i+z)-iy 

Tabs {Z, Zo) = Tabs ^ + ^) ~ ^ 



(A22) 

(A23) 
(A24) 



We have assumed that the photons are absorbed only by the Lya wedge, so the above 
equations are valid only for hu > 13.6 eV. To derive the photon fiux at energy lower than 
13.6 eV we have to take into account another hydrogen absorption lines. The radiation with 
energy below 13.6 eV start to be important only at lower redshifts (e.g. for rUs — 5 keV at 
z ~ 5) so it is not relevant in our calculations. 



B. Heating and ionization due to the sterile neutrino decays 

The photons from the decays of the sterile neutrino are mainly absorbed by neutral 
hydrogen atoms leading to their ionization. The ionization rate due to these photons is 
enhance almost 100 times due to additional ionization by the secondary electrons which 
deposit almost 1/3 of its energy into ionization (depending on the ionization fraction of the 
medium). The energy of the absorbed photons partially goes into ionization and partially 
into heating and excitations. We have adopted the approximation by Shull & Steenberg 
(1985), for which ionization rate and heating due to the photons from the decays of the 
sterile neutrino are respectively equal to 



A{z) = 



hu 



+ 



.He 



hu 
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X C,(l 



A-nau w) —, dv 

H hu 

th 



(Bl) 



where 



r, (z) 



+ 



Anan {i^) [hv - hug) dv 



AX 



hu 



hu 



CTHe H = 7.83 X 10~'^cm 



1.66 



.He\ 2.05 
th 



V 



-0.66 



,He\ 3-05" 
th 



V 



(B2) 



(B3) 



is the cross section for He ionization, hvg^ = 24.6 cV is its energy threshold; Ci = 0.3908, 
Ui = 0.4092, bi = 1.7592, Ch = 0.9971, = 0.2663, bh = 1.3163 (ShuU & Steenberg 1985). 



C. Cooling/heating function 

The H2 coohng function {Ah_h^, KH2-H2 ^^id ^form) given in Table 2 has been com- 
puted assuming that the CMB temperature is much smaller than the gas temperature. This 
approximation is not valid in our case since in the WDM model with keV sterile neutrinos 
the gas temperature can be very close or even lower than the CMB temperature because of 
relatively high ionization fraction. The situation is very similar in CDM model, where the 
gas temperature strictly follows the CMB temperature at least at high redshifts. 

In such conditions, the level population of H2 is strongly affected by stimulated emission 
and absorption. Molecules become an effective heating source for the gas, because the rate 
of coUisional de-excitation of the vibrational levels is faster than their radiative decay. In 
order to correct the coohng/heating function for this effect, we have to take the net cooling 
rate 

A,(r,7;) = A,(r)-A,(r,), (ci) 

where Aj represents either Ah-h2: ^H2-H2: or ^form- We also note that A^ {T,Tr) ^ Aj (T) 
for T > Tr. 

One should do the same for the Aj^_e (see Table 2). However, the coohng from the 
coUisional excitation of hydrogen atoms H and its subsequent spontaneous de-excitation is 
only important for the gas temperatures T > 10^. In our case the CMB temperature is well 
below 10^ K, so we can neglect the heating term for this cooling/heating mechanism. 
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For completeness we have to take also into account the fact that electrons can gain 
energy from the CMB radiation when T < 7]., so we are taking as the bremsstrahlung 
cooling/heating term 

In our calculations we have used the following cooling/heating function: 

AiT,Tr,z) = AH{T,Tr)+AHAT,Tr)+AcompiT,Tr) 

+ Ahrem{T,Tr)+Asiz). (C3) 
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Fig. 1. — The parameter range favored by the pulsar kicks, the current hmits from x-ray ob- 
servations (see text for discussion), and the benchmark points chosen for model calculations. 
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Fig. 2. — Evolution of temperature with redshift for different models: rris = 25 keV and 
sin^e = 3 X 10-12 (WDMl), = 15 keV and sin^ ^ = 3 x lO^^^ (WDM2), m, = 3.3 keV 
and sin^ 9 = 3 x 10^^ (WDM3). M is the cloud mass and z^r is the redshift of virialization. 




Fig. 3. — The same as in Fig. 2, but for the ionization fraction. 
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Fig. 4. — The same as in Fig. 2, but for the molecular hydrogen fraction. 
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Table 1. Reaction rates 



Reaction 



Rate (cm^ s ^ or s ^) 



Reference 



Ho formation via H channel 



e- + H -* H- +'y 
if + if- ^ i?2 + e- 
H- + J ^ H + e- 



ii^ 



H 



H ^ Ho 



'-2 -r JJ ^ -122 
H+ + 7 ^ ii + ii 



fei = 1.4 X 10-18T0-928 exp (- leioo) 

fe2 = 4.0 X io-9r-o-i^ 



fe3 = l.lxlO-iT2'i3exp(-8|P) 

H2 formation via channel 



logfc4 = -19.38 - 1.523 log T+ 1.118 (log T)2 
-0.1269 (log T)-^ 
fcs = 6.4 X 10""' 
fce = 2.0 X lO^T^^-SSexp 



(GaUi & Palla 1998) 
(GaUi & Palla 1998) 
(GalU & Palla 1998) 



(Gain & Palla 1998) 
(Galh & Palla 1998) 
(GalU & Palla 1998) 



CoUisional dissociation of H2 



H2 + if -» 3fi 



H2 + /i2 ^ 2ii + H2 



H2 + e- 



2ii + e- 



fc7,i, = 1.12 X 10-1° exp (-^^ 

=6.5 X IQ-'^T-i/^exp (^-^j [l -exp(- — 

logncr = 4.0- 0.416 log T4 - 0.327 (log r4)^ 

ks,L = 1.18 X 10-10 exp (^- 6.95X10-' ^ 
k8,H = 1.3 X 10-9 exp (-MSxloiJ 



1 + 



"«2 



logncr = 4.845 - 1.3 log T4 + 1.62 (logT4)^ 
kg = 1.3 X lO-l^T^exp (-^) 



(Shapiro & Kang 1987) 

(Palla, Salpeter & Stabler 1983) 

(Shapiro & Kang 1987) 
(Shapiro & Kang 1987) 



(Lenzuni, ChernoflF & Salpeter 1991) 



H2 + 7 ^ m 



H2 + 7 ^ 2fi" 



Photo-dissociation of H2 

2ii feio = 1.1 X 10* j (u) 

j (P) is the CMB radiation flux in erg cm— ^ s— 1 
at hp = 12.87 cV 
fell = /J? 47raii (u) j^dv 
14.159 eV, 

, 1 is the CMB intensity 



<^llH = TT^^riM + Ti^^fl M 



(T^j and <T°j axe the crossections for para- and ortho-H2 
respectively, y R! 3 is the ortho-, to paj:a-H2 ratio 



(Abel et al. 1997) 
(Abel et al. 1997) 



Ionization and recombination of H 



if + 7 ii'+ • 



fei2 = 47rCTi2 {v) ^du 
hug = 1*3.6 eV 



(Rybicki & Lightman 1979) 



ii+ + e- ^ if + 7 
H + e- -* H+ + 2e- 



ki3 = 8.4 X 10-11T-1/2T3-0-2 (1 + -1 
fci4 = 5.85 X 10-"Ti/2 (l + T^/^ye^p{-^ 



(GaUi & Palla 1998) 
(Haiman, Thoul & Loeb 1996) 
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Table 1 — Continued 



Reaction Rate (cm^ s ^ or s ^) Reference 



Ionization of H by the photons from the sterile neutrino decays 
(Lya ionization wedge and ionization by the secondary electrons) 

+ 7 -> i?+ + e- see text eq. (Bl) 



Note. — The binding energy of H2 is XH^i/k = 5.197 x 10* K. 
T and Tr = 2.725 (1 + z) are the gas and the CMB temperatures, 

respectively. is the CMB intensity. Ii^ (z) is the intensity of the 
photons from the sterile neutrino decays at redshift z. T„ = T/10" K. 
Xe (z), Zrec and Eq S'l'e the ionization fraction at redshift z, redshift 
of recombination and the energy of the photon produced in the sterile 
neutrino decays. 
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Table 2. Cooling/heating processes 



Cooling/heating function (erg cm ^ s ^) 



Reference 



H cooling 

Ah (T, Tr) = An-e (T) + Area (T) + (T) - TpD (Tr) 



(Haiman, Thoul & Loeb 1996) 
(Haiman, Thoul & Loeb 1996) 
(Haiman, Thoul & Loeb 1996) 



a) excitation of H by collisions with e 

An-e (T) = 7.5 X 10"" (l + T^/^) exp (-ilf« ) UeUH 

b) recombination to H 
Arec (T) = 1.04 X 10-^^ kiarienjj 

c) ionization of H by collisions with e~ 
Koniz (T) = 2.17 X 10-iifci4nen/f 

d) photo-ionization of H 

TpD (Tr) = inriH J^^ ^12 {i^) t ' ^""^h) 
/ii/ff = 13.6 eV 

H2 cooling 

{T,Tr)+A 

form 

{T,Tr) + (T) — Tpoi (Tr) — rpu2 (Tr) 

a) collisional excitation of H2 

i) excitation by collisions with H Ah-h^ {T, Tr) = Ah-h^ (T) - Ah-h^ (Tr) 

Ah-h, (T) = ScI'^^oha^ITt'e) ""H.nH (Galli & Palla 1998) 

ii) excitation by collisions with H2 ^H2-H2 {T,Tr) = AH2-H2 (T) - AH2-H2 i^r) 



^H2 — H2 



(T): 



A"^(n^O)+Ar^(,LTE) a"^ („^o)+A^'2 (iTB) / "^2 



A« (LTE) = A^ (LTE) + Af (LTE) 

log A" (n ^ 0) = -103 + 97.59 logT - 48.05 (logr)^ + 10.8 (logT)^ - 0.9032 (logT)" 



1+0.12T| 



-A — exp 



0.13 ' 



+ 3 X io-''4^^p|^_ih5i 



A^'-^2 tLTE) = — ^ 

^ ' "H,H 

Af ^LTE) = ^ [6.7 X 10-19 exp (-4^)" + 1.6 x 10-i« exp (-^)] 
A^2 ^ 0) = 0.25 [572^2 ^j^j + o.75 [fif ^ exp [- SJ^] 

Af2 (n 0) = 1.4 x 10-i2rV2 exp (- jJ^) exp ' 



7;^2 ^ 3 X 10-12+6.6 x IQ-i^Ts) |o 
|BJi/fe = BJo/it = 512 K, E^^/k = 5860 K 



T+1200 J 

276 J2 exp 



^10 \ c^v 
'~kT) ^10 

(sis, 



(HoUenbach & McKee 1979) 

(HoUenbach & McKee 1979) 
(Galli & Palla 1998) 

(HoUenbach & McKcc 1979) 

(HoUenbach & McKcc 1979) 

(HoUenbach & McKee 1979) 

(HoUenbach & McKee 1979) 

(HoUenbach & McKee 1979) 

(HoUenbach & McKee 1979) 



-A/orm (^) ^r) — ^form (T) A j- g^m (^r) 



3.53 



b) H2 formation cooling 
A/or-m (T) = 1.6022 X lo- 
ne. = 106r-V2/ {i.6x«exp [- + IAXH2 exp ["^^j } 

c) collisional dissociation of H2 

^diss (T) = 7.17145 X 10-12 (n^fcy + na^fes + nefcg) nH2 



Ife + 1.83^j^f|^) [1 - (1 + ^) nnn. (Shapiro & Kang 1987) 



fc2 4 




400 




T 





(HoUenbach & McKee 1979) 
(Shapiro & Kang 1987) 
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Table 2 — Continued 



Cooling/heating function (erg cm s ^) Reference 



d) photo-dissociation of H2 

^PDi (Tr) = 6.4 X W-^^kionH2 (Abel et al. 1997) 

ii) Ha + 7 ^ 

rpD2 {Tr) = 47rn«, [ ^ /J? (z.) ^ (hz. - huf^) di. + ^ a?, (z.) ^ (ftz. - ftz.°J di.] (Abel et al. 1997) 
hu^f^ = 14.159 eV, fet-"^ = 14.675 eV,y = 3 is the ortho-Ha to para-Ha ratio 

Compton cooling 

Acomp {T, Tr) = 1.017 X lO-s^T* (T - Tr) Ue (Haiman, Thoul & Loeb 1996) 



bremsstrahlung 

Afcrem {T,Tr) = Af^rem {T) ~ Af^rem {Tr) 

Air-em {T) = 1-42 X IQ-'^'^ Qf fT^/^ nenu+ (Haiman, Thoul & Loeb 1996) 

gff = 1.1 + 0.34 exp [- (5.5 - logT)^ /s] 

heating due to the photons from the sterile neutrino decays 



As (2) = -r„ {z) 



see text eq. (B2) 



Note. — T and Tr = 2.725 (1 + z) are the gas and CMB temperature respectively. Ji, is the CMB intensity. Iiy {z) is the intensity 

of the photons from the sterile neutrino decays at rcdshift z. T„ = T/10" K. ai (i/) and ki are given in the table 1. Xe (z), Zrec and 
Eq are the ionization fraction at redshift 2, redshift of recombination and the energy of the photon produced in the sterile neutrino 
decays. 



